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Abstract. The inclusion of running coupling effects in the BK-JIMWLK evolution 
equations considerably reduces the energy dependence of the saturation scale with 
' (— | ' respect to previous estimates based on fits to HERA data. We discuss how such 

Q-f slowdown affects the expectations for particle multiplicities in Pb-Pb collisions at 

O ■ LHC energies. Our prediction is based on the use of fc t -factorization and on the use of 

D | unintegrated gluon distributions taken from the numerical solutions of the BK equation 

i' with running coupling. We obtain a central value dN cfl j dr\{yfs = 5.5 TeV)\ v =o ~ 1390. 



> 

(S| ; PACS numbers: 21.65 Qr, 12.38 Mh 

The experimental results from RHIC [1] strongly suggest that heavy ion collisions 
at high energies probe QCD in the non-linear regime characterized by strong coherent 
fields and gluon saturation, the Color Glass Condensate (CGC) [2]. Thus, many bulk 
features of multiparticle production in RHIC collisions, such as the energy, rapidity and 
centrality dependence of multiparticle production, are succesfully described by models 
based in CGC physics [4,5]. With collision energies of up to 5.5 TeV, the upcoming 
^ ■ program in lead-lead collisions at the CERN Large Hadron Collider (LHC) is expected to 
provide confirmation for the tentative conclusions reached at RHIC and to discriminate 
between the different physical mechanisms proposed to explain particle production in 
high energy nuclear reactions (a review of predictions for heavy ion collisions at the 
LHC based on alternative approaches can be found in [3]). 

The phenomenological models in [4, 5] rely on the assumption that the saturation 
scale Q s a that governs the onset of non-linear effects in the wave function of the 
colliding nuclei is perturbatively large ~ 1 GeV at the highest RHIC energies. Next, 
gluon production is calculated via the convolution of the nuclear unintegrated gluon 
distributions (ugd's) according to ^-factorization [6] . Under the additional assumption 
of local parton-hadron duality, the pseudorapidity density of charged particles produced 
in a nucleus-nucleus collisions can be written as follows 

dN ch AnN c fd 2 p t f 2 / \h+Pt\ \ ( \kt-Pt\ 



dyd*b N*-lJ p2 y -""^.v^r^, 2 j r y» 2 

where tp denote the ugd's of projectile and target respectively, p t and y are the 
transverse momentum and rapidity of the produced particle, x 1>2 = (p t /y/s) e ±y , Q = 



2 



0.5 max {\p t ± k t |} and b the impact parameter of the collision. The lack of impact 
parameter integration in this calculation and the gluon to charged hadron ratio are 
accounted for by the constant C, which sets the normalization. 

The second basic ingredient of CGC models is the one of saturation of the nuclear 
udg's entering Eq. (TjQ) Using a relatively simple ansatz for the nuclear ugd's, and 
for symmetric collisions, the midrapidity multiplicity rises proportional to the nuclear 
saturation scale, which is assumed to grow as a power of the collision energy, ^/s: 

oc Q 2 sA ~ const ■ y/s X . (2) 



dyd 2 b y =o 

So far, the energy dependence of the saturation scale has been adjusted to the empirical 
value extracted from fits to DIS HERA data of saturation-based models [7], which 
yields A ~ 0.288. This has been largely motivated by the inability of the leading-log 
BK-JIMWLK evolution equations to reproduce experimental data. Here discuss how the 
recent advances in the calculation of running coupling corrections to the BK-JIMWLK 
evolution equations [8] allow to compute the energy evolution of the ugd's from first 
principles, getting a good agreement with experimental data. 

We start by solving the non-linear small- x evolution equation for the dipole-nucleus 
scattering matrix, S(Y,r), including running coupling corrections [8-10]: 

^^l = n[s\-s[S], (3) 

where r is the dipole size and Y = \n(xo/x). Explicit expressions for the kernel terms 
in Eq. ([3]) as well as a detailed explanation of the numerical method used to solve 
Eq. ([3]) are given in [9]. The initial conditions for the evolution are taken from the 
McLerran-Venugopalan model [11]. We classify them according to their initial value 
of the saturation scale, Qq. Importantly, the saturation scale extracted from solutions 
of Eq. ([3]) grows with energy significantly more slowly than the value extracted from 
empirical parametrizations of DIS data, as shown in Fig. [TJ On the other hand, the 
evolution speed A = din Q 2 s /dY turns out to be a function of Y, rather than a constant. 

We now use the exact solutions of Eq. (J3]) to calculate particle multiplicities 
according to Eq. (CQ). Thus, the nuclear udg's are now given in terms of the dipole 
scattering matrix evolved according to Eq. ([3]): 

<p(X,k) = j ^exp{tr-k}(l-S(Y,r)), (4) 

The relation between the evolution variable in Eq. ([3]) and Feynman-x of the produced 
particle is taken to be Y = ln(0. 05/0:1,2) + AF e „. Since the relevant values of Bjorken-x 
probed at mid-rapidities and ^s NN = 130 GeV at RHIC are estimated to be ~ 0.1-^0.01, 
the free parameter AY ev controls the extent of evolution undergone by the nuclear gluon 
densities resulting of Eq. ([3]) prior to comparison with RHIC data. Similar to [4] , large- 
x effects have been modelled by replacing (p(x,k) —>■ (p(x,k)(l — x) 4 . The running of 
the strong coupling in Eq. (DO), evaluated according to the one loop QCD expression, is 
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Figure 1. Left: Saturation scale as a function of rapidity corresponding to solutions 
of Eq. ((3]) (solid line) and setting A = 0.288 (dashed line). Right: Evolution speed, 
A = corresponding to the plot in the left. In both cases the initial saturation 

scale is Qo = 1 GeV. 



regularized in the infrared by freezing it to a constant value ctf r = 1 at small momenta. 
Finally, in order to compare Eq. with experimental data it is necessary to correct 
the difference between rapidity, y, and the experimentally measured pseudo-rapidity, 77. 
This is achieved by introducing an average hadron mass, m. Remarkably, the optimal 
value found in comparison with data, m ~ 0.25 GeV is in good quantitative agreement 
with the hadrochemical composition of particle production at RHIC. 

With this set up we find a remarkably good agreement with the pseudo-rapidity 
densities of charged particles measured in — 6% central Au+Au collisions at collision 
energies -Jsnn = 130 and 200 GeV, as shown in Fig. [2j The comparison with data [1] 
constrains the free parameters of the calculation to the ranges: Qo~0.75 -j- 1.25 GeV, 
m~0.25 GeV and 3> AY ev >0.5. Finally, the normalization constant is C ~ 0(1) in all 
cases. The smallness of AY ev indicates that the nuclear udg's probed at RHIC are in the 
pre- asymptotic regime [12]. With all the free parameters now constrained by comparison 
to RHIC data, the extrapolation to central Pb-Pb collisions at LHC energies, \/s = 5.5 
TeV per nucleon, is straightforward and completely driven by the non-linear small-x 
dynamics. We get 

dN Pb-Pb 

— c ± (y/J^= 5.5 TeV, 77 = 0) ~ 1290^ 1480, (5) 

ar] 

with a central value corresponding to the best fits to RHIC data ~ 1390. These values 
are significantly smaller than those of other saturation based calculations [4, 5, 13], 
~ 1700 -T- 2500, and compatible with the ones based on studies of the fragmentation 
region [14], which yield « 1000 1400. 
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Figure 2. Pseudo-rapidity density of charged particles produced in Au-Au 0-6% 
central collisions at y / s]vjv = 130 and 200 GeV and for Pb-Pb central collisions at 
V sjviv = 5.5 TeV. Data taken from [1]. The upper, central (solid lines) and lower limits 
of the theoretical uncertainty band correspond to (Qo — 1 GeV, Ay = 1), (Qo = 0.75 
GeV, AY" = 3) and (Q Q = 1.25 GeV, AY" = 0.5) respectively, with m = 0.25 GeV in all 
cases. 
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